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Abstract 

A new class of two- level non-regular fractional factorial designs is defined. We 
call this class an affinely full- dimensional factorial design, meaning that design 
points in the design of this class are not contained in any affine hyperplane in the 
vector space over F2. The property of the indicator function for this class is also 
clarified. A fractional factorial design in this class has a desirable property that 
parameters of the main effect model are simultaneously identifiable. We investigate 
^ , the property of this class from the viewpoint of D-optimality. In particular, for the 

VjD ! saturated designs, the Z)-optimal design is chosen from this class for the run sizes 

r = 5, 6, 7 (mod 8). 
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'. 1 Introduction 



In the literature on two-level fractional factorial designs, regular fractional factorial de- 
signs have been mainly studied both in theory and applications. The reason is that 
properly chosen regular fractional factorial designs have many desirable properties of be- 
ing balanced and orthogonal. In addition, the regular fractional factorial designs are easily 
constructed and used based on important concepts such as resolution and aberration. An 
elegant theory based on the linear algebra over F2 is well established for regular two-level 
fractional factorial designs. See p6] for example. 

On the other hand, non-regular designs have also been receiving attention of re- 
searchers over the years, in particular for some specific topics such as Plackett-Burman 
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designs, Hall's designs and mixed-level orthogonal arrays. However, it is very difficult 
to derive theoretical results for general non-regular fractional factorial designs. Suppose 
we have s controllable factors of two-levels and want to construct a fractional factorial 
design with r runs, where r is neither a 1/2-, 1/4-, 1/8-, . . . fraction of the 2* full factorial 
design. One approach to choose an r-runs design is to rely on various optimal criteria 
such as D-, A- or E'-optimality. However, the problem of characterizing these optimal 
designs for practically many values of (s,r) is combinatorially very difficult. See [28] and 
[7] for examples of D-optimal saturated designs. These works, in addition to enormous 
literature on Hadamard matrices such as j^, indicate the difficulty of the characteriza- 
tion of the D-optimal designs and the development of simple algorithms for obtaining 
them. For the other criteria such as A- or i?-optimality, see [33]. Another approach for 
the problem of optimal selection is various extension of the minimum aberration criterion 
to non-regular designs. See [13] and [39] for minimum G'2-aberration, [17] for generalized 
minimum aberration, ^6] for minimum moment aberration. 

In this paper, we give a new approach for investigating general non-regular fractional 
factorial designs from a theoretical viewpoint. We define a class of non-regular designs, 
which is derived naturally from the argument of identifiability of parameters. There are 
some other works considering such a classification of designs. [1] gives a classification 
of fractional factorial designs with simple structure in view of their indicator functions. 
The indicator function, ffist introduced in has become a powerful tool for studying 
general non-regular fractional factorial designs. See [S] and [2B] for example. Since our 
work is motivated by the indicator function approach, we also consider properties of the 
indicator function of the designs in our proposed class. 

The construction of this paper is as follows. In Section 2, we define a new class of the 
two-level non-regular fractional factorial designs and investigate its property. In Section 
3, we consider relations between the class and D-optimal designs. In particular we pay 
special attention to the saturated designs. Finally in Section 4, we give some discussion. 



2 Definition of an affinely full-dimensional factorial 
design 

Suppose there are s controllable factors of two levels. Let V be the 2* full factorial design 
with levels being —1 and 1. P is written as 



p = {-i,ir = {(xi,...,x,) \xi 



!}• 



A fractional factorial design JF (without replication) is a subset of V. Let r be the run size 
of J^. Therefore JF is a set of r points in V. Arranging the elements of J-" appropriately, 
we represent as an r x (s -|- 1) matrix M E { — 1, 
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where rriij is the level of the jth factor in the ith run. We call the left-most column of 
M the "0-th column" and write its elements as rriiQ = l,i = 1, . . . ,r. Note that M is 
the design matrix of JF for the main effect model. If we write jSj as the parameter of the 
main effect of the jth factor for j = 1, . . . , s, the linear model for the observation variable 
Y = (Yi, . . . , YrY is written as 

Y = M(3 + e, £~iV^(0,a2/), 
/3=(/3o,/3i,..., 

where Pq is the parameter for intercept. We put rriij = 1 for j = 1, . . . , s without loss of 
generality (by relabeling the two levels of each factor). We also assume that each column 
{mij, . . . , rrirj)', j = I, . . . , s, contains at least one —1. 

The indicator function / of is defined in [17J as a function on V such that 



/(x) 



1, ifxGJ^, 
0, ifxGl?\J^. 



Following [17], we define contrasts X/(x) = Yliei ^« P for / G V, where V is the set of 
all subsets of {1, . . . , s}. A fundamental fact in the theory of the experimental design is 
that {Xj, I G V} forms an orthogonal basis of the set of all real-valued functions on T>. 
The indicator function / of is then written as the polynomial form 

/(x) = 5^6,X,(x). (3) 

lev 

Since xf = 1 on X> for i = 1, . . . , s, the above expression is square-free and is unique on 
V. 

A regular 2*~^ fractional factorial design A gV is generated by k linearly independent 
generating relations 

X,,(x) = l,...,X,,(x) = l, (4) 

i.e. ^ = {x G I Xj^{x) = 1, i = l,...,k}. A contains 2*"^^ points. Note that the 
right hand sides of the relations (jlj) reflect the assumption rriij = 1 for j = 1, . . . , s. In 
general, we can take X/^(x) = —1 instead of X/^(x) = 1 in (jlj). This just depends on the 
labeling of two levels for each factor. From randomization viewpoint, given the labelings 
of two levels of each factor, it is desirable to choose 1 or —1 randomly for each generating 
relation in (jlj). 

Now we define a class of non-regular fractional factorial designs. 

Definition 2.1. A non-regular fractional factorial design T is called an affinely full- 
dimensional factorial design if there is no regular fractional factorial design A satisfying 
J-' G A. Conversely, a non-regular fractional factorial design JF is called a subset fractional 
factorial design if there is some regular fractional factorial design A satisfying T C- A. 

The above definition gives a new class of the non-regular fractional factorial designs. 
Any fractional factorial design with r > 2'^"^ is an affinely full-dimensional factorial 
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design. The merit of the definition might not be clear at a glance. One of the properties 
of the affinely full- dimensional factorial design is the simultaneous identifiability of the 
parameters in the model ([2]). Note that the least squares estimator of (3 is written as 
(3 = {M' M)~^M'y for the observation y. We call (3 is simultaneously identifiable if M'M 
is a non-singular matrix. We also note that, for regular fractional factorial designs, the 
singularity of M'M simply corresponds to a confounding relation of the main effects. On 
the other hand, for the non-regular designs, the relation is not obvious. We clarify this 
point in the following lemma. 

Lemma 2.1. // JF is an affinely full- dimensional factorial design, then M'M is non- 
singular, i.e., (3 is simultaneously identifiable. 



Proof. We show the contraposition. Suppose the columns of M are linearly dependent. 
This dependence relation is preserved by the following operation: subtract the left-most 
column (1, . . . , 1)' of M from the other s columns and then divide the s columns by —2. 
By this operation, 1 is mapped to and —1 is mapped to 1. Denote the resulting matrix 
by M = {mij}. Then = (1 — mjj)/2 for j > 1 and fhiQ = 1 for i = 1, . . . ,r. From 
the assumption of linearly dependence there exists some integer vector c = (cq, ci, . . . , Cs)' 
satisfying ^ 

Mc = 0. (5) 

Considering the modulo 2 reduction of ([5]), we have 

Mc = (mod 2), (6) 

where the odd elements and the even elements of c are replaced by 1 and 0, respectively. 
Here we can assume that there exists an odd element of c in since if every element of c 
is even then we can divide by the power of 2 in the factorization of the greatest common 
divisor of the elements of c. Moreover cq = 0, since the first row of M is (1, 0, . . . , 0). Then 
implies that there are even I's in {"fhij | = 1, j > 1} for i = 1, . . . , r, or equivalently 
there are even —I's in {rriij \ cj = l,j > 1} for i = 1, . . . ,r. Therefore X/(x) = 1 holds 
ioT I = {j \ Cj = lJ>l}. Q.E.D. 



From this lemma, we also have the following corollary. 

Corollary 2.1. Let M be a design matrix of a non-regular fractional factorial design. If 
M'M is singular, then the design is a subset fractional factorial design. 

The proof of Lemma 12.11 clarifies the geometrical meaning of the class. By the corre- 
spondence (1, —1) (0, 1), Xj{xi, . . . , Xs) = 1, (xi, . . . , Xs) eV if and only if 

CiXi -!-■■■ + CgXs = (mod 2), Xj = (1 — Xj)/2, j = 1, . . . , s, 

where 

' I 0, ifj^J. 
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However, as remarked earlier, design points can be chosen by Xj{xi, . . . , Xg) = — 1. There- 
fore, in general there may be a constant term Cq = or 1: 

Co + ciXi H h CgXs = (mod 2), xj = (1 - Xj)/2, j = 1, . . . , s. (7) 

In this sense, JF is a proper subset of no regular fractional factorial design if and only if 
the points of JF, considered as a subset of F2 = {0, 1}*, are not contained in any affine 
hyperplane of the form ([7]), i.e., we can form a basis of = {0, 1}* as the differences of 
the vectors in JF. This is the reason we call the class affinely full- dimensional. 

The question of determining whether a given design is an affinely full-dimensional is 
immediately read off from its indicator function. 

Lemma 2.2. Let T he a fractional factorial design and ^ be its indicator function. Then 
T is an affinely full- dimensional factorial design if and only if \bj\ < h^/^ for all / G P. 

Proof. Obvious from Proposition 4.2 and Corollary 4.3 of [17]. Q.E.D. 

3 D-optimality of the affinely full-dimensional facto- 
rial designs 

In Section 2, we define a new class of the non-regular fractional factorial designs, namely, 
affinely full- dimensional factorial designs. The affinely full-dimensional factorial designs 
have a desirable property that all the parameters are always simultaneously identifiable in 
the saturated model. Next problem of interest is whether this class includes good designs 
or not in view of various optimality criteria. In this paper, we consider D-optimality of 
the designs. For given r and s, the D-optimal design is the matrix M G { — 1, 
which maximizes | det(M'M)|. The D-optimal designs minimize the generalized variance 

3.1 D-optimal designs for the saturated cases 

First we consider the saturated cases, i.e., the cases of r = s + 1. In this case, the 
maximization of | det(M'M)| reduces to the maximization of | detM|, M G { — 1, l}'"^''. 
This problem is known as the Hadamard maximal determinant problem. Despite a century 
of works by mathematicians, this problem remains unanswered in general. 

To investigate the relation between the maximal determinant problem and the affinely 
full-dimensional factorial designs, we show a basic theorem. 

Theorem 3.1. A design with a design matrix M is affinely full- dimensional factorial if 
and only if det M is not divisible by 2^. 

We define M and M as in the proof of Lemma 12. 1[ The columns of these matrices are 
numberol from to s = r — 1 and the rows are numbered from 1 to r. Note that the first 
row of M is (1, 0, . . . , 0). This matrix has the following property. 
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Lemma 3.1. det M is an odd integer if and only if M is a design matrix of an affinely 
full- dimensional factorial design. 

Proof of Lemma 13.11 If M is a design matrix of a subset fractional factorial design or 
a regular fractional factorial design, there exists 1 < ji < ■ ■ ■ < jp < s, p > 2 satisfying 

p 

n ^^^^ = 1 
i=i 

for i = 1, . . . , r. From the definition of M, this is equivalent to 

p 

fhij^ = (mod 2) 

for z = 1, . . . , r. It means that M is a singular matrix in F2, i.e., det M = (mod 2). 

Q.E.D. 

Proof of Theorem 13.11 From the definition, we construct M from M as follows: 

• subtracting the left-most column of M from the other columns, 

rriij ^ iriij - niiQ, i = l,...,r, j = l,...,s, 

• dividing the columns except for the left-most column by —2, 

rriij ^ rriij /{-2), z = l,...,r, j = 1, . . . , s. 

Note that, after the first operation, no column of M is (0, . . . , 0)' from our assumption. 
Therefore (—2)^"^ det M = detM holds. From Lemma [3.11 we have proved the theorem. 

Q.E.D. 

Theorem 13.11 shows that whether a given non- regular fractional factorial design is affinely 
full- dimensional or not is judged from its determinant. Using this characteristic, we 
investigate the D-optimal designs. 

There are many literature reporting the solution of the Hadamard maximal deter- 
minant problem for specific r. The most basic result was given by Hadamard (^lOj) as 
det M < r''/^, where the bound is achieved only for Hadamard matrices. Paley has con- 
jectured that a Hadamard matrix exists for every r = (mod 4). The lowest order for 
which a Hadamard matrix is not yet obtained is r = 668 ([35]). For a design where the 
design matrix M is a Hadamard matrix, it is a subset fractional factorial design since the 
bound reduces to {AkY^^ = In fact, the product of the elements in each row is 

1 from the property of Hadamard matrix. Therefore we consider the case that r is not 
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Table 1: Summary whether the maximal determinant matrices are affinely full- 
dimensional factorial design (Yes) or subset fractional factorial design (No) for r = 
4, ...,99. The values of the maximal determinant are given as divided by 2^^^. The 
values that achieve the bounds ([2], [H], [15] and [15]) are underlined. 
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divisible by 4. According to the literature reporting the maximal determinant matrices, 
we investigate whether the D-optimal design matrix is affinely full-dimensional factorial 
or not for r < 100. The results are given in Table [TJ In Table [H "?" means that the 
maximal determinant matrix is not obtained (or the maximality is not proved) at present. 
We also show the values of the maximal determinant divided by 2*""^ for r = 1, 2, 3 (mod 
4). We underline the values, when the values achieve the known bounds such as [2] for 
r = 1 (mod 4), [13] and [35] for r = 2 (mod 4) and [15] for r = 3 (mod 4). We omit the ref- 
erence for the Hadamard matrices. For the other maximal determinant matrices, see [27] 
for r = 5; [33] for r = 6, 7; [IS] for r = 9, 10, 11; [3S] for r = 9, 10, 15; [?5] for r = 10, 26; 
[Tg for r = 11; [33] for r = 13, 25, 50; [12] for r = 13; [13] for r = 14, 18, 26, 30, 38; [ID] 
for r = 14,18; ^ for r = 15,74,82,90,98; [H] and [29] for r = 15; [31], [28] and [21] 
for r = 17; [50] for r = 18,26,30,38; [7] and [30] for r = 21; [4J for r = 25,61; [23] 
for r = 26,30,38,42,46,50,54,66; [3], [33], [3Z] and [31] for r = 41; [38] and [22] for 
r = 42, 66; [5l] for r = 50, 62; p] for r = 54; [ST] and [32] for r = 61; pj for r = 66; [5] 
for r = 86; [23] for r = 90. 

Table [H shows an interesting periodicity. We present the following conjecture. 

Conjecture 3.1. For r = 5,6,7 (mod 8), D-optimal design is affinely full- dimensional 
factorial. For r = 0, 1, 2, 3, 4 (mod 8 ), D-optimal design is subset fractional factorial. 

Though the authors do not succeed in proving this conjecture, we do not find any 
counterexample at present. We consider special cases instead. For the case r = 1 (mod 
4), ^ gives a bound 

det M < (2r - if'^r - if^^^^ (8) 

See also [T^ and [35]. This bound is achievable only if 2r — 1 is a perfect square. In 
fact, for r = 5,13,25,41,61 in Tabled! the maximum determinant achieves the bound, 
whereas the maximal determinant matrix has not been founded at present for r = 85. 
For the cases that the bound is attained, it is easy to prove the conjecture. 

Proposition 3.1. For r = 1 (mod 8) and the hound ^ is attained, the maximal deter- 
minant matrix is chosen as a design matrix for the subset fractional factorial design. On 
the other hand, forr = 5 (mod 8) and the bound ^ is attained, the maximal determinant 
matrix is chosen as a design matrix for the affinely full- dimensional factorial design. 

Proof. Write 2r — 1 = and r = 4/c + 1. Then it holds 

(2r-l)i/2(^_l)(-i)/2^^2^-ip/=_ 

Since m is an odd integer, the maximal determinant is divisible by 2*" if and only if k is 
an even integer, i.e., r = 1 (mod 8). Q.E.D. 

Similarly, for the case that r = 2 (mod 4), [^ and [3S] give a bound 

detM<2(r-l)(r-2)("-2)/2. (9) 
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In Table [H all the maximum determinants achieve the bound for r = 6, 10, . . . , 98 except 
for unsolved cases, r = 22, 34, 58, 70, 78, 94. For the cases that the bound is achieved, it 
is also easy to prove the conjecture. 

Proposition 3.2. For r = 2 (mod 8) and the bound ^ is attained, the maximal deter- 
minant matrix is chosen as a design matrix for the subset fractional factorial design. On 
the other hand, forr = 6 (mod 8) and the bound ^) is attained, the maximal determinant 
matrix is chosen as a design matrix for the affinely full- dimensional factorial design. 

Proof. For r = 8k -\- 2, k > 1, the bound, 

2(r - l)(r - 2)(^'^2)/2 ^ _ i^k^k^r+ik-i^ 

is divisible by 2"^ since A; > 1. For r = 8A; + 6, A; > 0, the bound, 

2(r - l)(r - 2)('^-2)/2 _ _ i)(2A; + l)^''-^^r-\ 
is not divisible by 2*" since both r — 1 and 2k -\- 1 are odd integers. Q.E.D. 

For the cases that r = 3 (mod 4), though a similar bound is given by [15], the maximal 
determinant matrix attaining the bound is not found at present except for r = 3. 

3.2 D-optimal designs for 4 and 5 factors cases 

Next we consider the non-saturated cases of s = 4, 5. We enumerate all the design matrices 
of the runs r < 10 and obtain the D-optimal designs. We investigate whether the D- 
optimal designs are the affinely full-dimensional factorial design or the subset fractional 
factorial design. The result is summarized in Table [2] and [31 Note that the factorial design 
is always affinely full-dimensional for r > 8, s = 4 since r > 2'^"^. 

We have not yet derived any clear relation between the D-optimality and the affine 
full-dimensionality like Conjecture 13.11 at present from the results in Table |2] and [31 



4 Discussion 

We classify non-regular designs by whether it is a proper subset of some regular fractional 
factorial design (namely, subset fractional factorial design) or not (namely, affinely full- 
dimensional factorial design). We also give a geometrical interpretation of each class, i.e., 
the affinely full-dimensional factorial design is characterized as the design with the points 
which are not contained in any affine hyperplane. One justification of our definition is 
the simultaneous identifiability of the parameters for the main effect model, which is one 
of the most important concepts in the theory of designed experiments. Therefore, if the 
main purpose of the data analysis is the estimation of the parameters for the main effect 
model, the strategy of considering the affinely full-dimensional factorial designs is useful. 
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Table 2: .D-optimal designs for 4 factors 



D-optimal 
design 



class 
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max I det M'M\ 
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"1 i i r 
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-1 1 1-1 
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D-optimal 
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regular FF 
X1X2X3 = 1 

215 
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Xi X2 X3 X4 

~ 1 1 r 

111-1 
11-11 
1-111 
1 -1 -1 -1 
-1111 
-1 1 -1 -1 
-1 -1 1 -1 
-1 -1 -1 1 

afRnely full-dim. 
none 
2^2 . 13 
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Xi X2 X3 X4 

~ 1 1 r 

111-1 
11-11 
1-111 
1 -1 -1 -1 
-1111 
-1 1 1-1 
-1 1 -1 -1 
-1 -1 1 -1 
-1 -1 -1 1 

afRnely full-dim. 
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2^2 . 3 • 7 
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Table 3: D-optimal designs for 5 factors 
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-1 


-1 


1 


-1 


1 


-1 


1 


-1 


1 


-1 


-1 


-1 


-1 




-1 -1 


1 


1 


-1 


-1 


1 


-1 


-1 


1 


-1 


1 


-1 


1 


1 




-1 -1 


1 


-1 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


-1 




-1 












-1 


-1 


1 


-1 


1 


-1 
-1 


-1 
-1 


1 
1 


1 

-1 


-1 
1 


class 


regular FF 






subset FF 






affinely full-dim. 




relation 


2:1X2X3 


= X1X4X5 


= 1 




X1X2XS = 


: 1 








none 






max|detM'M| 




218 










2I6 . 7 








2i4 . 72 
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It is an interesting topic to investigate whether various "good" designs are the affinely 
full-dimensional factorial designs or not. In this paper, we consider this problem from the 

concept of D-optimality. In particular, from the investigation of the D-optimal designs 
for the saturated cases, interesting mod 8 periodicity is suggested. Though the conjec- 
ture is not proved at present, it holds for all the proved £)-optimal designs under 100 
runs. We beheve that the conjecture interests researchers studying Hadamard maximum 
determinant problem. 

It is also important to consider other optimal criteria such as A- or i?-optimality. As 
an initial investigation, we have confirmed that the D-optimal designs of saturated designs 
for s = 4, 5 in Table 2 and 3 are also 74-optimal and £^-optimal designs. 
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